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Transport through semiconductor nanostructures is a quantum-coherent process. 
This review focuses on systems in which the electron's dynamics is balhstic and the 
transport is dominated by the scattering from structure boundaries. Opposite to 
the well-known case of the nuclear reactions, the potentials defining semiconductor 
structures are nonspherically symmetric and the asymptotic motion of the electrons 
is determined by the different potential levels in the contacts. For this special type 
of potential the mathematical foundations for the scattering theoretical description 
of the transport phenomena are presented. The transport properties of the system 
are then derived from the scattering matrix using the Landauer-Biittiker formalism. 
A rigorous analysis of the analytical properties of the S matrix leads to the most 
general resonant line shape described by a Fano function with a complex asymme- 
try parameter. On this basis the resonant and nonresonant contributions to the 
conductance and capacitance of the system are identified. 

I. INTRODUCTION 

A remarkable feature of the physics of semiconductor materials is the possibility of de- 
signing and manufacturing artificial structures in which the electrons can be confined in zero, 
one, and two dimensions^. The region of confinement, usually called quantum system, is 
coupled to contacts through tunneling barriers which assure the dominance of the quantum 
effects in the transport phenomena. Due to the coupling to the electrons in the Fermi sea of 
each contact the discrete electronic states of the isolated quantum system become resonant 
states^. They are special states of the continuum spectrum associated with a maximum 
of nonzero width of the electronic probability distribution density. The resonance width 
reflects the open character of the quantum system in semiconductor nanostructures and 
gives a measure of the coupling strength between the quantum system and contacts. The 
resonances can be directly seen in the transport properties of nanostructures giving reason 
to call the transport through this special type of mesoscopic structures as resonant-. 

In 1970 Easki and Tsu'^ have proposed for the first time a semiconductor structure where 
electronic transport proceeds via a resonant tunneling process and after that in 1973 they 
have also proposed a model^ to calculate the current density through the considered system. 
These were the first steps which have opened an extremely rich field for basic and applied 
research. The resonant transport has been investigated in a multitude of different mesoscopic 
semiconductor systemsii^: two dimensional electron gas, quantum point contacts, quantum 
wires and quantum dots. Many phenomena such as universal conductance fluctuations^. 
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the Aharonov-Bohm oscillations^, the quantum Hall effect^, the quantized conductance in 
ballistic point contactsiS, Coulomb blockade oscillationgAi, chaotic dynamics in quantum 
dotai^, and Kondo effect in single electron transistors-^^,, have been observed and discussed 
in the well-known theory pioneered by Landauer and Biittike r^'^i^^ . The widely successful 
application of this formalism shows that electron transport through a mesoscopic system is 
quite similar to scattering in nuclear or atomic physics^^. The most important difference 
is that for transport the asymptotic motion of the electrons is determined by the different 
potential levels in the contacts and therefore is not free as in the nuclear reaction theory. 
In turn, the spherical symmetry of the problem is broken and the methods to solve the 
scattering problem cannot be directly imported from the theory of nuclear reactions. The 
peculiarities of the scattering potential require a new theoretical description of the scattering 
phenomena appropriate for the transport. In this theory not only the values of the potential 
in the contacts should be taken into consideration, but also the strength of the coupling 
between the quantum system and the contacts. Recent research on Fano resonances in the 
conductance of a single electron transistor—, on fluctuations of the local density of states 
in the emitter^^ and on Luttinger liquid behavior in ballistic transport through quantum 
wires^^, has established the importance of the interaction of the quantum system with the 
contacts. As an effect of this interaction the line shape of the resonant profile is asymmetric 
and cannot be described by the common Wigner-Breit distribution anymore. 

The main aim of this paper is to present the mathematical foundations for the scattering 
theoretical description of coherent transport phenomena and to describe the relevant reso- 
nances for the transport properties. The paper is structured as follows: Sec. ITTl presents the 
method to determine the scattering matrix and the scattering functions associated to the 
scattering potential of the nanostructure. In the second part of this section the electronic 
charge and current densities are deduced in the frame of the Landauer-Biittiker formalism 
using the second quantization technique. Sec. IIIII presents an analytical theory of the quasi- 
isolated transport resonances, and the signature of these resonances in the conductance and 
capacitance measurements performed on semiconductor nanostructures is identified. The 
Fano functions with a complex asymmetry parameter arise as the most general resonant line 
shape and we provide explicit expressions for the parameters of the Fano profile. 

II. SCATTERING APPROACH TO QUANTUM TRANSPORT 

A. THE CONSIDERED SYSTEM TYPE 

The wave function which describes the electronic state of energy E is a solution of the 
Schrodinger equation 

^{E,f)=0, (1) 

where P is the momentum operator, Veff{r) is the effective potential energy which is a 
sum of the heteroj unction conduction band discontinuities, the electrostatic potential due 
to the ionized donors and acceptors, the self-consistent Hartree and exchange potentials due 
to free carriers and external potentials. The electronic states in mesoscopic systems are 
easily described within the effective mass approximation whose validity requires that the 
envelope function \E'(i?,r) be slowly varying over dimensions comparable to the unit cell of 
the crystal^S. 
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FIG. 1: A system relevant for transport properties: active region with size 2d sandwiched between 
source and drain contacts. The dimensions of the system are large enough to be extended to infinity 
{d < L^,Ly,Lz). 

We restrict our analysis to the case of semiconductor nanostructures presented schemati- 
cally in Fig. ^ in which the current flows only in one direction, from the source- to the drain 
contact, called transport direction; any kind of residual current is neglected. This means 
that the electrons are either bound in the lateral directions or freely moving in a plane per- 
pendicular to the transport direction, so that the mean value of the lateral current is zero. 
Consequently, the potential energy in which the electrons move should have a perpendicular 
component V±{r), which we assume independent of the coordinate in the transport direc- 
tion, i.e. V±{f) = V±{f±). For the structures which contain only a 2DEG and no lateral 
modulation V±{f±) = 0, whereas for quantum dots V±{f±) is considered either a parabolic 
or a square inflnite well potential energy. We assume, also, that except for this overall 
component V±{r±), there is no other lateral potential which can supplementary conflne the 
electrons in a certain region of the structure. Thus, the potential energy is separable. 



This special form of Veff{r) seems to be a strong restriction and it is not easy to imagine 
its validity for structures with complex geometry, such as a single electron transistor—, 
but it is the price payed for a good analytical description of the transport phenomena in 
semiconductor nanostructures. 

Due to the separable form of the potential, a particular solution of the Schrodinger 
equation ((T)) associated to the energy E = E± + e is 



where indexes the energy levels E'^ associated with the motion in the lateral directions. 
Here s is the degeneracy index for the electron motion in the transport direction. The 
functions (pi,{r±) and the energies E'^ are the solutions of the eigenvalue problem 
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The Hamilton operator of this problem is a Hermitian one and, consequently, the eigenfunc- 
tion set forms a basis for which the orthogonality and completeness relations 

Y.<P:{r^)Url) = S{f^-rl), (5) 

are valid. The index v is generic; for a free electron system in a plane perpendicular to 
the transport direction this index should be replaced by a 2D wave vector k^. The energy 
spectrum becomes continuous [E^ = fi'^k\/2m*). Thus, in Eqs. ^ 5uv' is transformed 
into 5{k± — k'j_) and the sum over u becomes an integral over k± (Ylu ~^ ^/(^^r)^/ dk±, 
where A is the lateral area of the system). Every energy level E± is degenerate with an 
infinite degeneracy. In the opposite limit, for a complete confinement of the system in the 
lateral directions, the energy levels are discrete and nondegenerate, and z/ is substituted by 
a quantum number. 

The z-dependent function in Eq. Q is a solution to the one dimensional eigenvalue 
problem 

^ :Pf + V{z) - e 
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^W(e,;.) = 0. (6) 



For nanostructures, the potential energy V{z) has some common properties: There is an 
active region which is small at the scale of the whole system. Inside this region, the electrons 
are elastically scattered from interfaces between different layers or between allowed and not 
allowed (depletion) domains in the case of a structured 2DEG. This region is also called 
scattering region and it is embedded between two contacts, which are practically semi- 
infinite homogeneous semiconductors. The electrons inside the contact regions move in a 
potential energy given only by the bulk conduction band structure. A relation between 
the surfaces of the scattering region and the interfaces between heterostructure layers or 
domains of the 2DEG can not be established. Each interface can cause elastic scattering 
and should be included in the scattering region. Thus, the planes z = ±d are chosen 
inside the homogeneous materials of the contacts, far enough from any interface and, of 
course, the choice is not unique. From the mathematical point of view, the potential energy 
V{z) is a function which varies considerably over small distances inside the scattering region 
and is constant outside, i.e. V{z < —d) = V\ and V{z > d) = V2, where ±d are the 
limits of the scattering region. Further we introduce the notations: Vm = min(Vi, V2) and 
Vm = max (1/1,1/2). 



B. S-MATRIX AND SCATTERING STATES 

The above considerations about the potential energy experienced by the electrons allow 
us to reduce the Schrodinger equation ((T]) to the effectively one dimensional eigenvalue 
problem 0. This problem is essentially a scattering problem for a particle with mass m*, 
in a potential with nonspherical symmetry, and it is usually solved in the scattering matrix 
approacbi^. 

In the stationary description of the scattering process it is assumed that the wave function 
of the system at large distances from the scattering region is a superposition of an incident- 
and a scattered wave^i. In our case the electrons are elastically scattered inside the domain 
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[— d, d\ and for all points outside this region the 2;-dependent part of the wave function is 
given by 

C V'*"(e, -d) exp[zA;i(z + d)] + V°"*(e, -d) exp[-iA;i(z + d)], z < -d 
i^{e,z)={ (7) 
[ ^*"(e, d) exp[-zA;2(z - rf)] + ^°"*(e, d) exp[zA;2(-2 - rf)], z>d 

where '?/'*"'/°"*(e, id) are complex coefficients and 



/c,(e) = ^^(e-K), s = l,2. (8) 

If e > then is a nonzero positive number according to the definition of the real 
square root function. However, if e < we have to take the first branch of the complex 
square root function, so that kg = i\ks\. At this stage we would like to point out the main 
difference between a spherically symmetric scattering problem like in the nuclear physics 
and a transport scattering problem: in the former only one wave vector for every value of 
the total energy is defined, while in the latter there are two wave vectors, fc^, s = 1,2, for 
each energy associated with the electron motion in the transport direction. If the materials 
of the source and drain contacts are different or when a current fiows through the structure, 
Vi and V2 are different and consequently ki ^k^. 

For a fixed energy e there are at most two independent solutions of Eq. (jH)) and, conse- 
quently, at most two expansion coefficients 

can be independent. We choose 
them (e, itd). The coefficients corresponding to the outgoing waves in Eq. ((Tj) should be 
expressed in terms of (e, ic?), 

7/'°"*(e, ) ~ ^^^^ \ #"(e, ) • 

The matrix S(e), defined for every value of the energy e associated with the motion of the 
electron in the transport direction, contains all the information about the collisions in the 
system^^ and is called scattering matrix. As a consequence of the existence of two wave 
vectors for each value of the energy e the scattering processes are described by a 2 x 2 
matrix even in the simplest case. 

For energies smaller than the absolute minimum of the potential energy, the solution 
given by Eq. ((Tj) should be zero and ^^"/°"*(e,±rf) = 0. If the energies are greater than 
min[l^(z)], the eigenvalue spectrum of the one dimensional Schrodinger equation © has a 
discrete part for \4i > e > min[V^(2;)] associated with bound states and a continuous part 
for e > Kn corresponding to the scattering states. We are interested in describing elastic 
processes which contribute to the transport phenomena and in this case the bound states 
do not play any role; Further we assume Vi^z^ > Vm- The continuous part of the energy 
spectrum has two important regions: e > Vm, for which two independent solutions of Eq. 
(jni) exist, and Vm < e < Vm with only one solution. For the degenerate energy levels there is 
not an unique way of defining the eigenf unctions and we further prefer to use the so called 
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scattering function^, 

ai^ T/^ f exp [iA;i(2; + (i)] + Sii(e) exp + rf)], z < -d 

^''\e.z) = %M\ , (10) 



2^ \S2iie)exp[ik2iz-d)], z>d 



9{e - V2) 



Suie) exp[—iki{z + d)], z < —d 



^(^)(e,z) = ^^^i^ . (11 



'^'^ I exp [—ik2{z — d)] + 822(6) exp [ik2{z — d)], z > d 

z) corresponds to a particle which comes into the scattering system from the left 
reservoir (s = 1) and is either transmitted or reflected. We can identify the element Sii(e) 
of the S matrix with the reflection amplitude, r^^\e)^ and the element S2i(e) with the 
transmission amplitude, t^^\e). ip^'^\e, z) corresponds to a particle coming from the right 
side of the system (s = 2). The reflection and transmission amplitudes for this particle are 
r'^^)(e) = S22(e) and ^'•^-'(e) = Si2(e), respectively. In the case e < Vs, the time translation 
invariance requires that |r'^'^)(e)| = 1, s = 1,2. 
The step functions 

^(e-K) = (i' s = l,2, (12) 

ensure that for \4i < e < Vm there is only one scattering function for each energy e. 

Generally, the wave functions are solutions of a differential equation and are defined up to 
a constant factor. In the above expression, we fix this constant, so that, in the limit of free 
particles [V{z) — > 0, d — »• 0], the scattering functions [Eqs. pOHlH) ] become eigenfunctions 
of the momentum operator P^. 

The matrix elements of S, i.e. the transmission and reflection coefficients, are formally 
determined by the continuity conditions of the scattering functions and their first derivatives 
at z = ±d, 

*(e) = -^e(e) [1 + S(6)] (13) 



^27r 

with the matrices (*(e)),,, = ^("')(e; (-l)'c/) and Bss' = 0{e - K)5ss', s,s' = 1,2, and 

*5(e) = --^^0(e)K(e) [1 - S{e)] (14) 



with (*s(e))^^, = (-l)** and K^^' = ks/{iT/2d)6ss', s,s' = 1,2, respectively. 

z={-l) = d 

We can use these conditions only after the calculation of the scattering functions inside the 
scattering region. 

In case of a general potential inside the scattering region, we can not represent the func- 
tions (e, z) by elementary functions between —d and d and it is necessary to expand them 
into a basis of eigenfunctions of a solvable problem. In the R matrix formalism the Wigner- 
Eisenbud basis is considered, i.e. the eigenbasis of the Hamilton operator corresponding to 
the closed counterpart of the studied scattering problem. 



Xiiz) = (15) 
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with the boundary conditions d xil dz\^^^^ = 0. Xi{^) a-^'s real functions defined only for 
z G [—d, d] and verify the orthogonality and the completeness relations, 



dz xi{z)xi'{z) = Sw, 



-d 



Y,Xi{z)xi{z') = 6{z-z'), 



(16) 
(17) 



1=1 



respectively. The R matrix method developed by Wigner and Eisenbud^^ and intensively 
used in the nuclear reaction theory2^, has become important in the last ten years also for 
describing transport phenomena24. 

We expand the wave functions inside the scattering region, z G [— c?, c?], in terms of the 
Wigner- Eisenbud functions 



ij^^\e,z) = J2<^'f\^)Xi{z) 



(18) 



1=1 



and the expansion coefficients a\'^^ are given by 



dz xi{z)ij^'\e,z) 



(19) 



for each s = 1, 2 and / > 1. To determine a^^\e) we use the Schrodinger equation © and Eq. 
(fT^ satisfied by the Wigner- Eisenbud functions. We multiply them by Xiiz) and ^^^\e,z), 
respectively, and the difference between the resulting equations is integrated over the interval 
[— (i, c?]. Performing an integration by parts in the kinetic energy term and identifying the 

integral on the right side with the coefficient a^i^\ defined by ()19|) . we obtain 



Xi{-d) 



2m* e — ei dz 



Xi{d) dtp'-'^ 



2m* e — El dz 



(20) 



Inserting Eq. ()2()j) into Eq. ()18j) we obtain the expression of the scattering functions between 
—d and as a function of their derivatives at the edges of the scattering area. 



i^^'\e,z) 



2d 

TT 



-R(e; —d, z] 



dz 



— R{e; d, z] 



s) 



dz 



where the R function is defined as 

R{e;z, z') 



2m* 2d 



1=1 



(21) 



(22) 



The above expression of the R function has the advantage of being dimensionless. Together 
with the continuity condition (fT^ the relation ()21|) gives the scattering functions inside the 
scattering region {—d<z<d)m terms of the S matrix. 



¥'\e,z) 



9{e - Vi) 2d 
V27V vr 
e{e - V2) 2d 



2n TT 



{tki [1 - Sn(e)] i?(e; -d, z) - zfc2S2i(e)i?(e; d, z)} , (23) 
{-^fciSi2(e)i?(e; -d, z) + [1 - 822(6)] i?(e; d,z)}. (24) 
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The values of the scattering functions on the surface of the scattering domain are then given 
by 



*(6) 

with the R matrix defined as 



e(e)R(e)K(e) [1 - S{e)] 



R{e; —d, —d) i?(e; d, —d) 
i?(e; —d, d) i?(e; d, d) 



(25) 



(26) 



The scattering functions are continuous at the edges of the scattering system and we 
exploit this condition to determine the relation between the R and the S matrix. Equating 
the two expressions (fT!^ and of the matrix ^ we infer 



2[l + 2R(e)K(e)]-' 



(27) 



for each value of the energy e in the domain of the scattering states. For Vm < e < Vm 
the above relation can only be used to calculate the matrix elements Sis and 82^ with s 
satisfying the condition e > V^, s = 1, 2. Using Eq. (jTTj) the elements of the scattering 
matrix can be effectively determined in terms of Wigner-Eisenbud functions and energies. 
Furthermore, the scattering functions, Eqs. pOllll|) and ()23H24p . can be evaluated in every 
point of the system, even inside the scattering system. 

As shown in the first part of this section, the elements of the S matrix defined by Eq. 
Q are in fact the reflection and transmission amplitudes for the particles coming from the 
left and from the right side of the system. However, for describing transport phenomena 
in semiconductors in the frame of the Landauer-Biittiker formalism^i, the reflection and 



transmission probabilities, R^^\t 



.(2) 



(e)| , T«(e) = W\e)\ and 



T(2)(, 



play a central role. We introduce further the current scattering matrijs24 
S(e) = Ki/2(e)S(e)K-i/2(e), (28) 



having the property that 



Sii(e) 



Si2(6: 



T(2)(, 



S2i(e) 



T(i)(e), and 



822(6) 



i?^^^(e). The relation ()27|) between the R and the 8 matrix can be rewritten into 



an equivalent form 

S(e) = l-2[l + iri(e)]"^ 
where the matrix of rank two is defined as 



00 , 



1=1 



with 



2m* 



s, s 



1, 2 



(29) 



(30) 



(31) 



for all / > 1. Per construction is a symmetrical matrix with real elements in the domain 
of the scattering states with e > Vm- This property leads to the unitarity of the current 
scattering matrix. 



e > V, 



Ml 



(32) 
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which includes the reciprocity relations, R^^\e) = R^'^\e) = R{e), T'^^^e) = T^'^\e) = T(e), 
and the well-known relation of the flux conservation, R{e) + T(e) = 1. 

The scattering functions ■ijj^^\e, z) form an orthogonal and complete system in the case 
of a potential energy V{z) which does not allow bound states^^. They are eigenfunctions 
of the Hamilton operator [Eq. (jH))] and the self-adjointness of this operator— allows us to 
consider the set as complete. 



POO 

J2 / degs{e){ij^^\e;z)yij'^^\e,z') = 6{z-z'), 



(33) 



where 

fji* 

defines the ID density of states. The orthogonality condition is written as 

dz (^W(6, z)y^P^''\e', z) = e{e - V,)5ss'S{e - e')/g,{e) (35) 
as demonstrated in Appendix El 

C. OBSERVABLES 

The conduction electrons in a semiconductor are considered as an electron gas embedded 
in a positively charged medium which maintains the overall charge neutrality of the system. 
The simplest approximation for the description of an electron gas is to neglect all interac- 
tions: the Coulomb interaction of the electrons with each other and the interaction of the 
electrons with the positive background. Every electron is then independent from other elec- 
trons and subject only to external forces. Each state of the system is described by the field 
operators \l/(r) and \I'^(r). Using the eigenbasis of the one-particle Hamiltonian associated 
with the Schrodinger equation (Q), we can represent ^(r) and \l/'^(r) in terms of creation 
and destruction operators for the states (z/, s, e) of the field, Cj,s(e) and c^^g^e), respectively. 
Thus, 

^(^1 = E / degs{e)Ur^)i^^'\e,z)cUe). (36) 
¥{r) = J2f degM€{r±) {¥'\e, z)y cUe), (37) 

with the ID density of states 5's(e), s = 1,2 defined by Eq. The operators ^(r) and 

'^\r) satisfy the fermion type anticommutation relations; consequently the creation and 
destruction operators have the following properties'^. 



{c..(e), c,.,.{e')} = 0, (38) 



4(e), clAe') ^ = 0. 
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The many-particle Hamiltonian of the electron system without mutual interaction, 

1 



H 



2m' 



-P' + V^{f±) + V{z) 



(39) 



can also be written in terms of creation and destruction operators, 



^ = 5Z/^^ (?.(e)^,(e)cL(e)c,.(e), 



(40) 



where E^{e) = E'^ + e is the energy of the single particle state (z/, s, e). An eigenvector |a) of 
the Hamiltonian (j4Up is completely determined by the occupation numbers nis^(e) of every 
single particle state (i^, s, e). 



\a) 



\n 



(«), 



(41) 



with clg{e)cus{€)\ni^\e)) = nl^\e)\nl^\e)) . According to the Pauli principle, at most one 
particle may occupy any fermion state so that these occupation numbers are restricted to 
the values and 1. 

The eigenvectors of the Hamiltonian, \a), form an orthonormal and complete system of 
vectors in the Hilbert space and each general state of the isolated electron gas can be given 
as a superposition of the pure states of H. According to the postulates of the statistical 
mechanics we suppose this superposition of pure states as incoherent^, so that the mean 
values of the electron charge and current densities in a mixed state are given by 



q{r) 



-2eTr 



p^^(r)^( 



and 



2eh^ 



Tr < p Im 



^^(f)V^(f) 



(42) 



(43) 



respectively, where p is the density matrix, e = |e| is the elementary charge and the factor 2 
comes from the spin degeneracy. In the definition of the electronic current density we have 
taken into account the conventional sign of the current. 

In the diagonal representation of the Hamiltonian the density matrix of the isolated 
system is also diagonal^^, paa' = Pa^aa', where Pa denotes here the probability of the system 
for being in the pure state |a) and satisfies the relation J2aPa = Tr[p] = 1. Inserting Eqs. 
(j36|) and (j37|) of the field operators in definitions (j42j) and (j43|) matrix elements such as 
{a\cl^{e)ciy/ gi {e')\a) occur. The special form of the orthogonality condition for the scattering 
states [Eq. yields 



Therefore, the electron charge and current densities become 

/>oo 



(44) 



(45) 
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and 



J[r) 



2eh X ^ 



de gsie) \ \(f)uif±)f Ira 



dz 



+ \^'^'\e,z)\ lm[<p:{f^)V,Mr±)]\nus{e), (46) 
respectively, where Uusi^) is the mean occupation number of the single particle state (z/, s, e), 



n,. 



Tr [p 



(47) 



and is the unity vector of the z axis. 

Further we analyze the component of the current density on the transport direction. 
Replacing in Eq. the scattering function by their expressions firU]) and (fTTj). and gs{e) 
by Eq. (jnil) we find 



Jz[r±,z) 



2e ^ 









2" 








/ de I n^i{e) 


1 - 


Sii(e) 






S2i(e) 


■1 


IVm I 















in the source contact region {z < —d) and 

. _ . 2e 



E 









2 






2" 


1 




nui{e) 


Si2(e) 




1 - 


S22(e) 






JVm 

















(48) 



(49) 



in the drain contact region {z > d). The nondegenerate energy levels between Vm and 
Vm do not contribute to the current because on the one side of the scattering region the 
corresponding scattering functions decay exponentially and on the other side the reflection 
coefficient is 1. For e > Vm the current scattering matrix S(e) is an unitary matrix and in 
turn the component of the current density in the transport direction can be written as 



Jz[r±,z) 



2e ^ 



deT{e)[n^i{e) - n,y2{^)] , \z\ > d 



(50) 



Vm 



where T(e) is the transmission probability. The current density is actually independent on 
z in the contact regions and has the same value in the source and drain contacts proving 
the current conservation through the system. 

To calculate the mean value of the occupation number ni,s{e) of the single particle state 
{i>, s, e) we need some considerations about the conduction electrons in nanostructures. The 
electron gas analyzed above extends in fact only in a domain whose dimensions are compara- 
ble with the phase coherence length of the nanostructure and are larger than the size of the 
scattering region. The rest of the heterostructure is usually a highly doped semiconductor 
or a metal and acts as a particle and energy reservoir for the considered electron gas. If the 
phase coherence length is large enough so that the electron gas can be taken as a macro- 
scopic system, the particle number and the energy of this system are practically constant in 
the limit of thermodynamic equilibrium. That means that the electron gas is quasi-isolated 
and the above expressions for the charge and current densities remain valid. The advantage 
of considering the system of electrons in contact with a particle- and energy bath is that the 
mean value of the occupation number is given by the Fermi-Dirac distribution function^^, 

nus{e) = fFDie + E'i-fi), (51) 

where /i is the chemical potential of the system at thermodynamic equilibrium, fixed by the 
doping in the contact regions of the heterostructure. The expression (j^T|) of the mean value 
of the occupation number assures that the current is zero at thermodynamic equilibrium. 
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D. LANDAUER-BUTTIKER FORMALISM 



To calculate the transport properties in the Landauer-Biittiker formalism^^, electrons 
can be thought of as two noninteracting Fermi-gases: First, the electrons coming from the 
source contact which occupy the single-particle scattering states with s = 1 according to the 
Fermi-Dirac distribution function fpD^E — fii), 



(52) 



where fii is the chemical potential of the source contact. Second, the electrons coming from 
the drain contact which occupy the single particle states indexed by s = 2 according to the 
Fermi-Dirac distribution function fpD^E — (12), 



nu2{e) = fpnie + E"^ - 7x2) 



(53) 



= A^i — GUsd being the chemical potential of the drain contact. The potential difference 
/ii — fi2 results from an externally applied drain-source voltage, Ugd- 

Using these basic assumptions of the Landauer-Biittiker formalism and the normalization 
condition (0) of the functions 0i/(rj_) we determine from Eq. the electronic charge 

Q = J dfq{r) inside the scattering system. 



Q 



2em* 



E 



Vs 



P (f) 



(54) 



where 



PJe) 



dz lip 



(55) 



is the particle probability distribution. For a free electron system in a plane perpendicular to 
the transport direction index u should be replaced by a 2D wave vector k± and the electron 
charge in the scattering region becomes 



Q = -A- 



TT 



m 
1? 



^de^ 
. Ue) 



dEfFD{E-fi,), 



(56) 



where E denotes the total energy of the electron and A is the lateral area of the system. It 
is the great advantage of using the R matrix formalism to enable the analytical calculation 
of -Ps(e) inside the scattering system, as shown in Appendix ^ Its expression is given by 
Eq. ()A7|) and using the expressions ()13II14|) for ^ and ^'s we obtain 



PJe) 



27r 



Tm 



for e> Vf 



^dS 

de 



+ 



(57) 



Pii^) 



2tt 



-Im 



9i{e) 



S* 
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,^^ dSn ^ Sn{e) 



de 



h{e) 



1 h{e) 
47r|A:2(e)|2 



521(e) 



(58) 
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for e < Vm in the case Vi < V2, and 



^2(6) 



2tt 



-Im 



92{e) 



n* ( ^ dS22 5*22(6) 



de 



hie) 



1 hje) 
An |A;i(e)|2 



Suie) 



(59) 



for e < Va/ in the case V2 < Vi. The above relations demonstrate that the electronic charge 
accumulated inside the scattering system is only given by the elements of the scattering 
matrix. 

To describe the transport properties of nanostructures we also need to analyze the current 
through the system, I = f dr^jz{r±, z) as a functions of the applied drain-source voltage. 
Using Eqs. (|52II53|) we obtain from Eq. (j^UI) the current as 



POO 

■ / de T(e) V [fM^ + EI- /ii) - ^(e + E^- /X2)] 

■JVm 



(60) 



for \z\ > d, where T(e) is the transmission probability characterizing the scattering region of 
the system. This is the current between the source and drain contacts. We can also define 
a current in the lateral directions as results from Eq. PU)) . but this current is zero for the 
systems considered in this paper. For the heterostructures with total lateral confinement, 
the energy levels E'^ are nondegenerate and the time translational invariance allows us to 
choose the function 4>u{r±) as real. It results immediately that Im [0* (ri)Vrj^0i/(rj_)] = 
and the lateral component of the current density vanishes. In the opposite limit, for a free 
electron gas in a plane perpendicular to the transport direction we can define a 2D wave 
vector kj_ and 4'u{^±) ^ 4'{k±,'r±) = e^^-^^^/27[. It follows in a straightforward manner that 
Im [0^(rl) Vr=*x0iy(rj_)] k±/ (27r)2 and the lateral component of the current density becomes 
zero. 

For determining the conductance it is convenient to write the expression (|6U|) of the 
current into an equivalent form 



2e 



POO 

/ dE [fME - fii) - fpDiE - fi2)]y2TiE - E''^)eiE ~ E-'^), 

JVm 



(61) 



where the integration is made over the total energy of the electron. This expression has the 
advantage that it directly yields the conductance. Each energy level i?^ defines a channel for 
the electron transport^'^ and it is usually said that the channels which satisfy the condition 
< E\_ < Ep are open channels because they contribute to the current. The 6'-function in 
Eq. (jUT|) serves to remove the channels with exponentially decaying wave functions in the 
contacts, the so called closed channels for transport. In the linear response regime {Usd 0, 
Vi = V2 = 0) and for low temperatures (T — > 0), we can expand the Fermi-Dirac function 
/i?£)(i? — /X2) in a Taylor series around E — fii, and thus obtain^^ from Eq. 



2e 



(62) 



where T(e) = Si2(e 

dimensional scattering problem Eq, 
in the source contact. 



is the transmission probability determined solely by the one- 
^ and Ep is the Fermi energy of the electron gas 
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III. RESONANCES IN TRANSPORT 

The notion of resonances, representing long-lived intermediate states of an open system 
to which bound states of its closed counterpart are converted due to couphng to continuum, 
is one of the most fundamental concepts in the domain of quantum scattering^. On a formal 
level resonances show up as poles of the scattering matrix occurring at complex energies 
(^0 = Co ~'^r/2, where eo and T are called position and width of the resonance, respectively^. 

The causality condition leads to analytic properties of the scattering matrix when the 
energy or the wave vector of the electron are extended to complex values'^^. To calculate 
the poles of the scattering matrix we exploit these properties which are briefly summarized 
in the following. S{e) is a meromorphic function on the two-sheeted Riemann surfaces - 
one for positive values of the imaginary part of the wave vector and the other for negative 
values - with a branch point at e = and a cut from to oo. Bound states poles lie on the 
negative real axis of the "physical sheet" and, except for them, S{e) is an analytical function 
on the physical sheet^. The resonance poles lie on the second, "unphysical sheet" and come 
from the possible zeros on the first sheelj^. As usually done, we assume these poles to be 
simple^. 



e -plane (physical sheet) £ -plane (unphysical sheet) 





cut 


cut 


X X— 

Bound state 
poles 




>^ X 

Resonance poles 



FIG. 2: Left side: "Physical sheet" of the complex energy plane corresponding to the values of the 
wave vector k = \J 2m*e/fi'^ with Im[A;] > 0. Poles of the scattering matrix associated with the 
bound states. Right side: "Unphysical sheet" of the complex energy plane corresponding to the 
values of the wave vector with Im[A;] < 0. Poles of S'(e) associated to resonances. 

Our development of the resonance theory starts with the R matrix representation of the 
S matrix which is particularly well suited to the description of narrow resonances^. As 
demonstrated in Sec. Ill Bl the two matrices are related to each other through Eq. (j29p . It 
is obvious that the resonance poles correspond to the zeros of the denominator function in 
the expression (j^^ of the S matrix, 

det[l + in(eo -ir/2)] = 0. (63) 

According to the definition pup the 17 matrix elements have single poles for each e = e/, 
/ > 1, i.e. for each eigenenergy of the Hamiltonian corresponding to the scattering system 
isolated from the contacts. The Wigner-Eisenbud energies ei are real and are associated to 
the bound states. If the system is coupled to contacts the bound states are transformed into 
scattering states (at most with a few exceptions depending on the particularities of V{z)) and 
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the poles should migrate in the lower part of the "unphysical sheet" of the complex energy 
plan. Further we demonstrate that the Wigner-Eisenbud energies, e/, are not solutions of 
Eq. ()63|) and therefore they are not singularities of the S matrix. The matrix O is split into 
a resonant part which contains (e — ex)~^ and a regular matrix fix, 

I = 1 
I ^ A 

This decomposition and the property of ujx to have the determinant zero lead to 

det[l + in] = l^l^^Mllvxie) (65) 
e - ca 

with ^A(e) = — ^Tr[a;A(l + ^^^a)^^] and Vx{e) = det [1 + if^A] • Now we can express the 
S- matrix in terms of fix and a; a, 

S(.) = (66) 

e - eA - £x{e) 

where the matrix Zx{e) is defined as Zx^e) = (e — eA) [—1 — det[f2] + z ($7 — fi^)] 
with ri~ = r2~^det[0]. In principle, the solutions of the equation Vx{e) = can not be 
associated with resonance poles of the S matrix. As follows from its definition Sx{^) becomes 
infinite for 'Dxie) 0. We cannot argue that T>x and det[l + ifl] vanish at the same points 
and it is convenient to include Vxi^e) in Zx{e). 

The matrix Za and the function Sx are related to each other through the unitarity re- 
quirement for the S matrix which gives 

ZxZ{ = Z{Zx=\e-ex-Sx{e)\\ (67) 

for e > Vm- 

The representation of the S matrix in Eq. (j66p is an exact reformulation of Eq. (j29p and 
has the advantage of directly yielding the equation 

eo-ex- Sxieo) = (68) 

to determine the positions eo = eo — «r/2 of the poles in the complex energy plane. It is 
obvious that Eq. has no singularities for real energies in the interval (eA-i, ca+i). The 
coupling to the contacts leads to a nonzero imaginary part of the resonance energy and the 
stronger the coupling the larger the difference between the Wigner-Eisenbud energy eA and 
the resonant energy cqa, A > 1. Each Wigner-Eisenbud energy was associated a resonance 
energy. 

In the following we analyze the narrow transport resonances for which F is a small quantity 
and the line shape is given by an asymmetric Fano function. Therefore, as a basic assumption 
for our theory of the resonant Fano line shape, we require the validity of the linearization 
of Sx and implicitly Za in a domain of the complex energy plane that includes the pole 
eo and the part of the real axis which contains the transmission peak associated with the 
resonance, i.e the resonance domain. To obtain the line shape of the resonance we employ 
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a formal expansion of the S-matrix as given in Eq. (j66|l in a Laurent series around the pole 
Co = eo — and neglect the derivatives up to the second order of £\ and Za at the points 
eo and ecr^^. Thus, the S matrix has the form 



1 Z,(eo)-.r/2^| 

S e 



e - eo + iT/2 i _ ds^ 



de 







de \ e= 




1 dSx 




^ de 






e=eo 



(69) 



Eq. ()66|1 ensures that the S matrix is an analytic function in the resonance domain excepting 
the pole, condition which is required for the existence of the Laurent series. After the 
linearization of Ex in Eq. we find 1 — d£/de\^^^ = (eo — ex — £x{eo)) Ar/2 and then 



the expression of the S matrix inside the resonance domain becomes 

:S(eo) — Sbg 



S(e)-z . + S,g, (70) 

e + 1 



where e = 2(e — eo)/r and 

iT/2 d!Ly 



eQ- ex - i^A(eo) de 



(71) 



e=eo 



Eq. (fTOj) has a standard for m^^i^^i^^i^^ , but we can provide here an explicit expression for the 
nonresonant component of the scattering matrix, S^^. The first term in Eq. (|70|l represents 
the resonant part of S. For each element of the matrix S it is seen from Eq. ()70|) that the 
resonant part undergoes a phase change of tt when the energy passes the resonance. In 
general, this produces a change between constructive and destructive superposition of the 
resonant and the nonresonant part. Therefore, an asymmetric line is obtained. 

The resulting expression (j7(J|) preserves the unitarity of the scattering matrix only in 
linear order e^^: ^ 

SS^ = S^S ~ 1 + ((5 - 1)-^. (72) 

We expect that our approximation is valid as long as the second term on the right hand side 
of Eq. ()72j) is small compared to 1, i.e. the deviation of S from unitarity is small. This way, 
for each maximum we can estimate the range of validity for our approximation through the 
requirement (<5 — e^/ (e^ + 1) -C 1, i,j = 1, 2. 



IV. APPLICATIONS 



A. CONDUCTANCE THROUGH A QUANTUM DOT 

As a first application of our resonance theory of transport we calculate the conductance 
through a quantum dot embedded in a quantum wire-^. The potential energy experienced 
by electrons has two components: the z-independent lateral confinement potential energy 
V±{f±), which provides the one-dimensional character of the structure, and V{z) a double 
barrier potential separating the quantum dot from the rest of the system (Fig. El)- 

In the linear response regime {Vgd — >■ 0, l^i = V2 = 0) and for low temperatures {T —>■ 0), 
the conductance through the dot is given by Eq. (jU^ as a superposition of transmission 
curves T[Ep — E\/), where v corresponds to an open channel for the transport. 
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In experiments G and therefore T is probed at different energies by varying the voltage of 
an additional plunger gate. In case of lateral tunneling, this additional gate is a top gatei^^ 
and, in the case of vertical tunneling, it is a side gate^^. We use the following idealization 
for the total potential in presence of a varying gate voltage: The external potential created 
by the charges at the gate is screened out completely in the heavily doped contacts {\z\ > d) 
so that the total potential and Ep remain unchanged. In the scattering area {\z\ < d) 
the total potential can be idealized for small variations of the gate voltage as a varying 
potential energy offset ef/g, so that V{z) — > V{z) — ellg. As shown in Fig.Elthe transmission 
probability of the double barrier system depends on the gate potential Ug and therefore the 
conductance G varies with U„. 
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FIG. 3: Middle: Total potential energy of a double-barrier structure, V{z) — ellg, for a value of the 
gate voltage which corresponds to a maximum in conductance, Ug = Uq = 14.28mV. The potential 
steps of height —eUo at \z\ = d {d = 16nm) result from the voltage applied to the plunger gate. 
In dotted lines the wave functions at the resonant energies ej(C/o), i = 1,2. The energies 

Ep — E'^ (solid lines) are drawn in the contact regions, where E'^, v = 1,6 correspond to an infinite 
2D quantum well with the dimensions 16nm x I6nm. Left: Positions of the resonance poles of the 
S matrix in the complex energy plane. Right: Transmission T{Uq; e) vs. energy calculated for the 
scattering potential energy V{z) — eUo: numerical calculation (solid lines) and Fano approximation 
(dashed lines). 

As illustrated in Fig. |21for small energies (e < max[V"(z)]) the transmission is generally 
small and may have some isolated peaks at ej. The peak position depends on the total 
potential energy experienced by electrons and therefore ej = eiiUg). The line shape of a 
transmission peak is deduced from Eq. ()70|) and is given by 

T(e) = 



(S(e)) 



12 



Tr 



bg- 



e2 + 1 



(73) 



18 



where Ti 



bg 



112 



2 



is the background transmission. The right hand side of Eq. ()73|) is a 
Fano distribution with a complex asymmetry parameter— 

, = .(sm)Js.,);;. (74) 

As shown in Fig. Elthe Fano approximation for the hne shape of the transmission gives a very 
good description even for the asymmetric peaks as long as the peaks are quasi-isolated. The 
background matrix Sbg in Eq. (fTOj) is often assumed to be absent, leading to a symmetrical 
Wigner-Breit line shape of the transmission^^, but it is obviously that such a description 
fails in the case of an asymmetric peak as the second peak in Fig. IHl 

If the Fermi energy is smaller than the maximum of the barriers and if the external 
voltages Ug are quite small, only the isolated peaks participate to the transport. According 
to Eq. ()62|) the conductance has a maximum each time the transmission of an open channel 
has a peak, 

e,,{Uo) = Ep ~ (75) 

Thus, to each conductance maximum at Ug = Uq a pair index {uq, iq) is assigned, where z/q 
is the channel index and io the number of the maximum in the curve T{Uq\ e). As shown 
in Fig. ini each index iq can be associated with a pole of the scattering matrix and we can 
conclude that a maximum in conductance is characterized by a resonant channel z/q and 
a resonant pole of the S matrix, On this basis one can split the contribution to the 
conductance into a coherent resonant part, 

Gc{Ug) = ^T{Ug;Ep-El°), (76) 
and a slowly varying noncoherent part 

GNciUg) = ^ E ^(^s' - ^^)' (77) 

in which the absolute squares of the transmission coefficients are added without phase- 
information. In the above expression only the contributions of the open channels are counted. 
Fig. m illustrates the conductance through the dot and its coherent and noncoherent contri- 
butions around a maximum at Ug = Uq. 

In case of a narrow conductance peak the R-matrix representation of the scattering matrix 
allows for obtaining an explicit dependence G = G{Ug) in the vicinity of the maximum, 
Ug = Uq + 6U. By the transformation V{z) — > V{z) — eUg the Wigner-Eisenbud energies e; 
[Eq. (jl5|) ] become e; — eUg but the Wigner-Eisenbud functions xi remain unchanged. Then, 
according to Eq. (jHUj) the matrix 0(f/c,; e) can be approached by f2(f/o; e + e6U). Thus, 
from Eq. follows that T{Ug] e) ~ T{Uq; e + e5U). Inserting condition (f73|) in the above 
relation we obtain the coherent and noncoherent contributions to the conductance around 
the maximum, 

Gc{Ug) ~ ^^T{Uo-^ e,, + e5U), (78) 

and 

G^ciUg) ^ ^ E ^(^0' + ^I" -EI + e6U), (79) 
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respectively, where e^^ is the position of the tQ-th. transmission maximum for the potential 
energy V{z) — cUq. 

Using Eq. ()73|) the coherent contribution to the conductance in the vicinity of the reso- 
nance is obtained as a Fano function with the complex asymmetry parameter q defined by 



G 



bg- 



[e + Re(g)]'+ [Im(g)]^ 

2;2 



(80) 



Here e = 2e{Ug — f/o)/rjg is a function of the plunger gate potential f/g, the resonance 
position Uq and the resonance width calculated as the imaginary part of the iQ-ih pole of 
the scattering matrix for the potential V{z) — cUq. The background coherent contribution 
is related to the background transmission through Ghg = {2e^ /h)Thg. 

As shown in Fig. |3] the noncoherent part of the conductance varies slowly inside the 
resonance domain and we can expand this function in a Taylor series for Ug ~ Uq. One 
obtains 



GNc{Ug 



2^ 
h 



J2 T{Uo; + El" 



(81) 



'El 



As illustrated in Fig. |1[ Eqs. (jHUjl and (jHT|) provide a very good description of the two 
contributions to the conductance of the dot. If the overlap of the peaks is very small, the 
noncoherent contribution to the conductance can be considered even as constant^^. 
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FIG. 4: Left: Conductance vs. gate voltage around the maximum at Uo = 14.28 mV: complete 
calculation Eq. (solid line) and approximative value G = Gc + Gnc given by Eqs. (|5n|) and 
H81|) (dashed line). Right upper part: Coherent contributions to the conductance around maximum: 
complete calculation Eq. H76() (solid line) and Fano approximation Eq. (jHUj) (dashed line). Right 
lower part: Noncoherent contributions to the conductance around maximum: complete calculation 
Eq. H77|) (solid line) and linear approximation Eq. H81j) (dashed line). 

We can conclude that the most general form of a conductance profile is a Fano line with 
a complex asymmetry parameter superposed on a linear function which usually increases 
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the asymmetry of the profile. The Fano function arises from the coherent superposition 
of contributions to the S matrix coming from different poles and the linear function is 
determinate by the nonresonant channels whose contributions add noncoherently. Therefore 
we have demonstrated that the asymmetry of the profile in conductance does not necessarily 
involve the coupling between two different channels as in the usual scenario to explain Fano 
resonances^. In addition, the asymmetry of the conductance profile is not necessary the 
asymmetry of the transmission peak. 

In our theory of the resonant transport all parameter necessary to describe the conduc- 
tance profile can be calculated microscopically and we can evaluate separately the coherent 
and the noncoherent parts of the conductance. The two contributions arise naturally in our 
formalism of the coherent transport through open systems. We do not exclude the existence 
of the incoherent processes in nanostructures which can contribute to the conductance, but 
we have fond that the coherent processes contribute coherently and noncoherently to the 
conductance. 



B. CAPACITANCE OF QUANTUM MIS-TYPE HETEROSTRUCTURES 



Further, we analyze a quantum MIS (metal - insulator - semiconductor)-type heterostruc- 
ture. The considered AlAs/GaAs structure"^^ consists of a sequence of layers grown on a 
GaAs bulk material given by, first, n — GaAs layer as a back contact, second, intrinsic 
GaAs layer as a spacer, third, a short period Alri-Gai^^As / GaAs superlattice as a blocking 
barrier and finally a metallization as a top gate. The band structure. Fig. resulting from 



drain 




FIG. 5: Band diagram from self-consistent calculation^MS for the MIS-type AlAs/GaAs het- 
erostructure. The structure parameters correspond to the experiments in Ref,^^. 

the self-consistent calculations in the Hartree approximatio n^°i^° shows for large enough pos- 
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itive gate voltages Ug a potential well at the interface between the GaAs spacer and the 
blocking barrier, where a field induced two-dimensional electron gas (2DEG) develops. 

The blocking barrier is now assumed to suppress charge transfer completely and we 
consider the limit of small frequencies. Then, the tunnel capacitor in Fig. becomes 
equivalent to a simple plate capacitor The charge corresponding to the left plate is 
distributed in the region —d < z < Zf, (see. Fig. ^ and can be evaluated using Gauss' law, 
Qi = —Sn{dV/dz){z = Zb). Here Zb is the interface between the barrier and the GaAs- 
spacer layer, S is the area of the sample, and k is the dielectric permittivity. The top gate 
of the structure acts as the right plate of the capacitor and the charge on it is Q2 = —Qi- 
The capacitance is then readily found as 



C 



dQi 




dQ2 


dUg 




dUg 



(82) 



In Fig. ini we compare the experimental C-V-curve^^ with the results of our model^S. 
For the numerical calculations, we use the parameters corresponding to the experiments in 
Ref.— . Because the work function of the metal contact is not precisely known we shift the 
theoretical voltage scale Ug with respect to the experimental one Uq {Ug = 0.7011^ + Ug), 
so that the centers of the steps coincide. Ug = Q corresponds to the fiat band configuration. 
It can be seen that for such a system the C-V curve takes the form of a broadened step 
located between a low voltage, Ug < U_, and a high- voltage, Ug > f/+, plateau. We define 
the gate voltage Uc at the center of the step by (PC /dUg = 0; U- and [/+ correspond to the 
gate voltages where \dC/dUg\ takes half of its maximum value. 
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FIG. 6: C-V curve, experimental data (filled circles^ and numerical calculation (solid line). 

To understand the variation of the capacitance with the gate voltage in the step domain 
we perform the pole analysis. For the scattering potential of the biased structure. Fig. [7| 
illustrates the pole energies with the real part in the energy domain of the occupied states 
and the scattering functions at the resonance energies corresponding to these poles. This 
figure shows that there are two types of poles: the first one is associated with the scattering 
states of the electrons confined in the region close to the back contact. The corresponding 
charge distribution varies very slowly with the applied biaa^^*^. Therefore, these resonances 
do not play any role in the capacitance variation with Ug. The second type of poles, with a 
smaller imaginary part, corresponds to the states of the electrons localized at the interface 
between the GaAs spacer layer and the blocking barrier. In our case there is a single pole 
with the resonance energy below the chemical potential /^i and for which \il)'^^\e, z)\'^ has 
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FIG. 7: Left: Position of the resonance poles in the complex energy plane determined by Eq. 
Right: Potential energy V{z) (solid line) for Ug = 16 mV and \ip^^\€, z)p at the resonance energies 
(dotted line). The upper plots present in detail an energy domain around f^i. Fmax marks the 
height of the wide and shallow barrier which is formed between the back contact and the blocking 
barrier, V^a.^ = max_d<z<zt,[y{z)]. 

a maximum in the region of the potential quantum well^. We denote the complex energy 
of this pole by eo = eo — «r/2. The imaginary part of this pole is smaller than that of the 
other poles and smaller than the gap between two adjacent resonance energies. The charge 
accumulated in the quantum well varies strongly with Ug and, consequently, the changes of 
the capacitance can be directly connected to the resonances of the second type. A quantity 
which gives information about the contributions of different poles to the charge accumulated 
in the system is the probability distribution Ps^e) [Eq. EH))]. In Fig. |Hl -Pi(e) is plotted for 
two values of the applied bias: Ug = 11 meV, corresponding approximatively to the center 
of the capacitance step and Ug = 25 meV, corresponding to the high voltage plateau of the 
capacitance curve (Fig. In the middle part of Fig. ISlthe resonance energy eo and the 
energies eo ± r/2 are presented as a function of Ug together with the energy of the maximum 
of -Pi(e), emax, and the energies e± where -Pi(e) takes half the maximum value. For comparison 
the chemical potential and the height of the wide and shallow barrier which is formed between 
the back contact and the blocking barrier, Vmax = iiiax_d<2<zjl^(2;)], are also plotted. In this 
situation Knax is below the chemical potential /ii and a classically allowed channel for the 
electrons is opened. For small positive values of Ug the maximum of -Pi(e) lies in this channel. 
With increasing Ug a potential quantum well is formed at the interface between the spacer 
layer and the blocking barrier and -Pi(e) has a sharp maximum which is located in an energy 
range with no classically allowed connection between the electrons in the quantum well and 
the back contact. So we can conclude that for positive values of Ug -Pi(e) has a pronounced 
maximum centered on eo with a well defined width (~ F). The charge accumulation in 
the field induced quantum well of the quantum MIS type semiconductor is characterized 
mainly by the resonance associated with the pole eo — ir/2. This resonance changes its 
character from an intermediate resonance to a quasibound state. For small applied biases 
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P,(e) [nm] U„ [mV] P,(e) [nm] 

FIG. 8: Left part: Probability distribution Pi(e) around the resonance energy eg for Ug = 11 mV. 
Right side: Probabihty distribution Pi{e) around the resonance energy eo for Ug = 25 mV. Middle 
part: The energy Cmax of the maximum of Pi{e) (solid line), and the energies e± at which Pi{e) 
takes half the maximum value (dotted lines), the resonant energies eo (filled circles) and the energies 
eo ± r/2 (circles). Shaded area: energies which correspond to the classically allowed channel 

(;Ui > e > Frnax)- 

{U- < Ug < U+) this resonance has the character of an intermediate resonance^: i) its 
energy lies in the classically allowed range (Fig. |H1 middle part) and, therefore, is a Fabry- 
Perot type resonance; ii) it is located in the space between the contact and the region where 
an isolated 2DEG is formed at large values of Vg (Fig. 3b.) Ref.~); iii) its shape is strongly 
asymmetric (Fig. |H1 left part). With increasing Ug the intermediate resonance turns into 
a quasibound state which, in contrast to the intermediate resonance, is connected to the 
back contact only via the tunneling effect. The resonance line of -Pi(e) narrows and tends 
to become symmetric (Fig. |S[ right part). The excellent quantitative agreement between 
the experimental data and the modeled capacitance curve. Fig. IHl demonstrates that the 
electronic states, which are important for the capacitance step and the subsequent high- 
voltage plateau are derived from a single resonance. 

V. CONCLUSIONS 

We first analyze the electronic scattering states in semiconductor nanostructures and, 
on this basis, the transport properties of the electron gas: conductance and capacitance. 
The considered systems are mesoscopic and the transport is supposed ballistic. The physics 
of these systems is often dominated by resonances which are broadened due to the open 
character of the systems. We develop a theory of resonant transport in semiconductor 
nanostructures valid for all coupling regimes between the quantum system and contacts, i.e. 
for open systems as well as for almost-closed quantum systems. To separate weakly from 
strongly energy dependent contributions on the scale of the width of the resonance we employ 
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a representation of the S matrix in terms of the R matrix. We find a simple procedure to 
determine the poles of the S matrix in the complex energy plane, which gives the position 
and the width of the resonances. Numerical results show the remarkable accuracy of our 
procedure which is applied to a great variety of resonances types: quasibound states, Fowler- 
Nordheim, and Fabry-Perot resonances. In case of separate resonances, the linearization of 
the weakly energy dependent part in S matrix yields an analytical expression of the line 
shape which is a Fano profile with a complex asymmetry parameter. The resonant theory 
of transport is applied to study the conductance through a quantum dot embedded in a 
quantum wire and the capacitance of a 2DEG formed at the interface between the spacer 
layer and the blocking barrier in a MIS-type semiconductor heterostructure. 



APPENDIX A: ORTHOGONALITY OF THE SCATTERING FUNCTIONS 

The main aim of this section is to demonstrate Eq. (j35j) using the expressions (jlUj) . 
flTT|) . and (fTH|) of the scattering functions and the relation between the S and the R 
matrix. For an effective calculation of the scalar product between ip^^^e; z) and ^/'^^'•'(e'; z) it 
is necessary to split the integral into two contributions, the first one given by the scattering 
domain and the second one given by the lateral regions (outside the scattering domain), 

dz [^(^)(e;^)]>(^')(e';z) = /-/(e, e') + /rf/(e, e'). (Al) 

Further we are interested in evaluating the above integral only in the cases in which the 
energies e and e' are in the same domain: e, e' > Vm or < e, e' < Vm- If the two 
energies belong to different domains, the corresponding eigenfunctions are orthogonal due 
to the self-adjointness of the Hamilton operator in Eq. ® and consequently Eq. (jH^jl is 
fulfilled. 

On integrating over the lateral regions such terms occur as 



where V denotes the Cauchy principal part. 



(A2) 



77 > 

and after a laborious calculation we obtain 

ItPie, e') = e{t - Vs)e{e' - Vs>) {5 {kl{t) - k,{t')) [SuSw + SL(e)Si,,(e')] 

+5{kl{e) + A;i(e')) [Su^iAe') + SL(e)5i.'] 
+5{kl{e) - k2{e')) [62s62s' + SL(e)S2.'(e')] 
+6 {k;{e) + k^ie')) [52sS2s'(e') + SL(e)52s']} /2 
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To evaluate the delta functions and their coefficients we have to discuss the cases e, e' > Vm 
and Vm < e, e' < Vm separately. For e, e' > Vm we obtain 9{e — Vs)0{e' — Vs') = 1, 
6{k*{e)-k,{e')) = 6{e-e') / g^ie) with g,{e) defined by Eq. s = 1,2 and 6{k*{e)+ks{e')) = 
0. Using the relation between the S- and the S matrix and the unitarity of S [Eq. ijH^ ] the 
first term in the expression ()A4|) becomes 6{e — e')5ss' / Qsi,^)- In the case < e, e' < Vm it 
is found that e{e - - K') = 5(fc;;(e) - /c^(e')) = 5(e - e')/c/^(e) with c/^(e) 

defined by Eq. (jMjl and (5(fcj^(e) + fcm(e')) = 0- Taking also into account that for considered 
energies Smm(e) = 1, we obtain the same result as in the first case. In the last two terms 
of Eq. (jSH) we substitute [^^(e)]^ - ^^(e') = (e - e')2m*/Ti^, s = 1, 2, which directly follows 
from the definition JHl) of kg and /f5*(e, e') becomes 



9s[e) 



(A5) 



where the matrix Mfe, e') is defined as 



M(e,e') 



2m* e 



for e 7^ e' and 



2m* 



(A6) 



(A7) 



for e = e'. 

To calculate Il^J^{e,e') we use the expression (fTS|) of the scattering functions inside the 
scattering region and find 



/.".'(f.t'i 



2rfV 



TT y 



Ur^(6)) (*5(e')),v / M(e;(-l)X^)i?(6';(-iyrf,^)(A8) 



The integral in the above expression of /^"/(e, e') can be easily calculated if the R functions 
are replaced by their definition ()22j) and the orthogonality condition of the Wigner-Eisenbud 
functions [Eq. (fTHj)] is used. After that a very simple trick, 



V- 



where the Cauchy principal part is defined by Eq. ()A3|) . allows us to write 



(A9) 



dzR{e; (-l)M, z)R{e'- (-l^rf, z) 



2m* 2d e — e 
Inserting Eq. ()A10|) in Eq. ()A8|) we obtain 



-d 

[R{e'; i-iyd, i-iyd) - R{e; (-l)X (-1^^)] • 



M,,,(e,e') 



(AlO) 



(All) 
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Using the expressions of the integrals inside and outside the scattering region [Eqs. ()A11|) 
and ()A5|) . respectively] the orthogonality condition [Eq. of the scattering functions 

follows from Eq. ()A1|) . 



^ Ferry D. K.and Goodnick S. M. 1999, Transport in Nanostructures, Cambridge University Press, 
Cambridge. 

^ Kukulin V. I., Krasnopolsky V. M., and Horacek J. 1989, Theory of Resonances - Principles 

and Applications, Kluwer Academic Publishers, Dordrecht. 
3 Esaki, L. and Tsu R. 1970, IBM J. Res. Develop. 14, 61. 
^ Tsu R., and Esaki L. 1973, Appl. Phys. Lett. 22, 562. 

^ Kelly M.J. 1995, Low-dimensional Semiconductors, Clarendon Press, Oxford. 

^ Ando T., Arakawa Y., Furuya K., Komiyama S., and Nakashima H.(Eds.), 1998,Mesoscopic 

Physics and Electronics, Springer Verlag, Berlin. 
^ Fowler A. B., Hartstein A., and Webb R. A. 1982, Phys. Rev. Lett. 48, 196; 
^ Biittiker M., Imry Y., Landauer R., and Pinhas S. 1985, Phys. Rev. B 31, 6207. 
9 Biittiker M. 1988, Phys. Rev. B 38, 9375. 

10 van Wees B. J. et al. 1988, Phys. Rev. Lett. 60, 848; Szafer A. and Stone A.D. 1989, Phys. Rev. 
Lett. 62, 300. 

11 Meirav U., Kastner M. A., and Wind S. J. 1990, Phys. Rev. Lett. 65, 771. 

12 Jalabert R. A., Stone A. D., and Alhassid Y. 1992, Phys. Rev. Lett. 68, 3468. 

1^ Goldhaber-Gordon D., Shtrikman H., Mahalu D., Abusch-Magder D., Meirav U., and Kastner 
M. A. 1998, Nature, 391, 156. 

14 Landauer R. 1957, IBM J. Res. Develop. 1, 223; Landauer R. 1987, Z. Phys. B 68, 217; Fisher 
D. S. and Lee P. A. 1981, Phys. Rev. B 23, 6851; Biittiker M. 1986, Phys. Rev. Lett. 57, 1761; 
Biittiker M. 1988, IBM J. Res. Develop. 32, 317; Stone A. D. and Szafer A. 1988, IBM J. Res. 
Develop. 32, 384. 

15 Biittiker M. 1992, Phys. Rev. B 46, 12485. 

16 Humblet J. and Rosenfeld L. 1961, Nucl. Phys. 26, 529. 

1'^ Gores J., Goldhaber-Gordon D., Heemeyer S. , and Kastner M.A., Shtrikman H., Mahalu D., 
and Meirav U. 2000, Phys. Rev. B 62, 2188. 

18 Schmidt T., Konig P., McCann E., Fal'ko V. I., and Haug R. J. 2001, Phys. Rev. Lett. 86 , 276. 

19 Lai S., Rao S., and Sen D. 2001, Rhys. Rev. Lett. 87, 026801. 

20 Bastard G., Brum J. A., and Ferreira R. 1991, Sol.State Phys. 44, 229. 

21 Taylor J. R. 1972, Scattering Theory: The Quantum Theory of Nonrelativistic Collisions, John 
Wiley & Sons, Inc., New York. 

22 Wigner E. Rand Eisenbud L. 1947, Phys. Rev. 72, 29. 

23 Lane A. M. and Thomas R. G. 1958, Rev. Mod. Phys. 30, 257. 

24 Smrcka L. 1990 , Superlattices Microstruct. 8, 221; Wulf U., Kucera J., Racec P. N., and 
Sigmund E. 1998, Phys. Rev. B 58, 16209; Alhassid Y. 2000, Rev. Mod. Phys. 72, 895. 

25 Racec E. R. 2002, Ph.D. Thesis, University of Technology Cottbus. 

26 Huang K. 1965, Statistical Mechanics, John Wiley & Sons, Inc., New York. 

2^ Biittiker M., Imry Y., Landauer R., and Pinhas S. 1985, Phys. Rev. B 31, 6207. 

28 Racec E. R. and Wulf U. 2001, Phys. Rev. B 64, 115318. 

29 Racec P. N., Wulf U. and Kucera J. 2000, Solid-State Electronics 44, 881. 



27 



Racec P. N., Racec E. R., and Wulf U. 2002, Phys. Rev. B 65, 193314. 
Fyodorov Y. V. and Sommers H. J. 1997, J. Math. Phys. 38, 1918. 

Bohm A. 1993, Quantum Mechanics, Springer, New York, Chap. 18 Resonance Phenomena, 
van Kampen N. G. 1953, Phys. Rev. 91, 1267. 

Nockel J. U. and Stone A. D. 1995, Phys.Rev. B 51, 17219; Nockel J. U. and Stone A. D. 1994, 

Phys. Rev. B 50 17415. 

Biittikcr M. 1988, Phys. Rev. B 38, 12724. 

Landau L. D. and Lifschitz E. M. 1977, Quantum Mechanics (Non-Relativistic Theory), Perg- 
amon, Oxford. 

Alhassid Y. and Attias H. 1996, Phys. Rev. B 54, 2696; Alhassid Y. and H. Attias H. 1996, 

Phys. Rev. Lett. 76, 1711. 

Kastner M. A. 1993, Phys. Today 46, 24; 

Tarucha, S., Austing D. G., Honda T., van der Hage R. J., and Kouwenhoven L. P. 1996, Phys. 
Rev. Lett. 77, 3613. 

Racec P. N., Racec E. R., and Wulf U. 2001, Computational Materials Science 21, 475. 
Racec P. N. 2002, Ph.D. Thesis, University of Technology Cottbus. 

Dolgopolov V. T., Shashkin A. A., Aristov A. V., Schmerek D., Drexler H., Hansen W., Kotthaus 
J. P., and Holland M. 1996, Phys. Low-Dim. Struct. 6, 1. 
Cavicchi R. E. and Silsbee R.H. 1988, Phys. Rev. B 37, 706. 



